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Abstract 

I present results for two-loop soft anomalous dimensions for heavy quark production 
which control soft-gluon resummation at next-to-next-to-leading-logarithm (NNLL) accu- 
racy. I derive an explicit expression for the exact result and study it numerically for top 
quark production via e + e~ — * tt, and I construct a surprisingly simple but very accurate 
approximation. I show that the two-loop soft anomalous dimensions with massive quarks 
display a simple proportionality relation to the one-loop result only in the limit of vanish- 
ing quark mass. I also discuss the extension of the calculation to single top and top pair 
production in hadron colliders. 



The calculation of higher-order corrections to heavy quark cross sections and differential 
distributions is needed for increasing the accuracy of the theoretical predictions. Given the 
importance of the top quark [I] to electroweak and Higgs physics, the top quark cross section 
at the Tevatron and the LHC, both in top pair [2] and single-top [5] production modes, as well 
as at a future e + e~ collider needs to be calculated with the highest possible precision. Bottom 
and charm quark cross sections are also important in understanding QCD and calculating 
backgrounds to new physics. 

Soft-gluon emission is an important contributor to higher-order corrections, particularly 
near partonic threshold. The soft-gluon corrections can be formally resummed to all orders in 
perturbation theory. The resummation arises from the factorization of the cross section into a 
hard-scattering function H and a soft function S that describes noncollinear soft-gluon emission 
in the process [I]. The evolution of the soft function is controlled by a process-dependent 
soft anomalous dimension Tg. The resummed cross section can be written in moment space 
as a ICS (N) = exp[Ei(N,/j,)}H(ij,)S(N)exp[f™ /N (dfi/iJ,)2T s (ij)], with fi the factorization scale, 
m the heavy quark mass, and N the moment variable which is conjugate to a kinematical 
variable measuring distance from partonic threshold. The exponents resum collinear and 
soft radiation from any massless partons in the scattering. The logarithms of N exponentiate. 
At leading-logarithm accuracy we determine the coefficients of the highest power of In N at each 
order in the strong coupling a s , at next-to-leading logarithm (NLL) accuracy we also determine 
the second highest power of In N, etc. 

For heavy quark pair, such as ti, hadroproduction the soft anomalous dimensions are ma- 
trices in color space due to the complicated color structure of the scattering processes, and 
were calculated at one loop for both qq — > ti and gg — > ti channels in Ref. [3], thus allowing 
resummation at next-to-leading logarithm. From the calculation in [4] one can also derive the 
corresponding result for heavy quark production at e + e~ machines, a simple function related 
to the cusp anomalous dimension of Ref. [5] . For single top production at hadron colliders the 
one-loop soft anomalous dimensions were presented for the t, s, and tW channels in Ref. [6]. 
To increase the accuracy of the resummation to next-to-next-to-leading logarithms (NNLL) for 
all these processes one needs to derive the soft anomalous dimensions at two loops. 

The calculations of soft anomalous dimensions involve diagrams with eikonal (Wilson) lines 
representing the heavy quarks. The eikonal approximation is valid for describing the emission 
of soft gluons from particles in the hard scattering and leads to a modified form of the Feynman 
rules and thus of diagram calculations (see e.g. [4-9]). When the gluon momentum goes to zero, 
the quark-gluon vertex reduces to g s T F v^/v ■ k, with g s the strong coupling, v a dimensionless 
velocity vector, k the gluon momentum, and Tp the generators of SU(3) in the fundamental 
representation. Below we calculate explicitly the soft anomalous dimension through two loops 
for the process e + e~ — > ti, study its properties in detail and then discuss the extension to 
processes with more complicated color structure, i.e. single top and top pair hadroproduction. 

The eikonal diagrams are calculated in dimensional regularization with n = 4 — e and in 
Feynman gauge in momentum space. The one- loop soft anomalous dimension, Vg , can be 
read off the coefficient of the ultraviolet (UV) pole of the one- loop diagrams in Fig. 1. Writing 
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Figure 1: One- loop diagrams with heavy-quark eikonal lines. 
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we find the one-loop expression 
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with C F = 4/3, (3 = ^1 - Am 2 /s, and s = (p e + + p e -) 2 . 

We now continue with the two-loop diagrams. A few results for some of these diagrams 
have appeared in [HI El [TOj . In Fig. 2 we show graphs with vertex corrections and in Fig. 3 
graphs with heavy-quark self-energy corrections. There are three more graphs, not shown, with 
vanishing contributions: two graphs with gluon loops involving a four-gluon vertex, and one 
graph involving a three-gluon vertex with all three gluons attaching to a single eikonal line. 
The calculations are challenging due to the heavy quark mass and involve multiple complicated 
integrals and delicate separations of infrared and ultraviolet poles (which by construction of 
the soft function are opposites of each other [1]). The analytical structure involves logarithms 
and polylogarithms. It is understood in the results below that for each diagram we include 
the appropriate one-loop counterterms for the divergent subdiagrams (but note that graph 
2(b) does not have a divergent subdiagram). Typically there are large cancellations of terms 
between the individual diagrams and their respective counterterms. 

We now present the UV poles of the kinematic terms for the two-loop diagrams, which later 
will have to be combined with color and symmetry factors, and an overall factor a 2 /n 2 . The 
sum of the UV poles of diagrams 2(a) and 2(b) is minus one-half the square of the UV pole of 
the one-loop diagram of Fig. 1(a): I 2a + I 2b = [(1 + (3 2 ) 2 /(8/? 2 )](-l/e 2 ) ln 2 [(l - f3)/(l + (3)}. 
The UV poles of the crossed diagram, Fig. 2(b), are 
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The diagram in Fig. 2(c) represents corrections from quark, gluon, and ghost loops. The 
quark-loop contribution is I 2cq = [n/(l + (3 2 )/(Q(3)][l/e 2 - 5/(6e)] ln[(l - (3)/(l + with n f 
the number of light quark flavors, while the gluon plus ghost contribution is I 2c g = (5/24) [(1 + 
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2 )/0][l/e 2 - 31/(30e)] ln[(l - 0)/{l + (3)}. Diagram 2(d) is given by 
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and the result for diagram 2(e) is its negative: he — ~hd- The diagram of Fig. 2(f) gives 
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Diagram 3(a) represents two graphs with overall Cp color factor, I, 
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and another graph / 3a2 = 1/e 2 — l/(2e). The quark- loop contribution to graph 3(b) is hbq = 
(n//3)[l/e 2 — 5/(6e)] while the gluon plus ghost loop contribution is hbg = (5/12) [1/e 2 — 
31/(30e)]. Note the relation d/d0[20I 2c / {I + 2 )] = [-2/(1 - 2 )}hb- Finally diagram 3(c) 
gives J 3c = [(1 + 2 )/{20)}{-l/e 2 ) ln[(l - 0)/{l + 0)}. 

Combining the above kinematic results with color and symmetry factors, the contribution 
of the diagrams in Figs. (2) and (3) to the two-loop soft anomalous dimension is 
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where we have cancelled out factors of a 2 /ir 2 on both sides of the equation. On the right-hand 
side of Eq. ([5]) in addition to Y K S , which appears in the coefficient of the 1/e pole, there also 
appear terms from the exponentiation of the one-loop result and the running of the coupling, 
with O = (11/3)Ca — 2rif/3, Ca = 3, which account for all the double poles of the graphs. 
From Eq. ([5]) we solve for the two-loop soft anomalous dimension: 



,(2) 



K Ca 
2 2 



3 

{1 + 2 ) 



+ C F C, 



(1+0 2 



20 



40 

1-0 " 
1 + 0, 

In 



C 



' 1-0 " 
Li 2 



+ In 



1-0" 



1 

3 

0" 



In 2 



1 + 
' {1-0? ^ 
,(l + 0) 2 , 
l-0\ 



(2 



+ -ln 2 

3 



0" 



[1 + 



2\2 



1 + 



In 



1 + 0) 

qi+0? " 

V 4/? , 



S0 2 
6 



1 + 0/ 
-Li 3 

1-0 " 
1 + 0, 



+ C2 
'{1- 



0? 



{1 + 0? 



+ C3 
-0) 



1 + 0? 



(6) 



(2) 

We have written the two-loop result in Eq. 



in the form of a term which is a multiple 



of the one-loop soft anomalous dimension plus additional terms. The well-known two-loop 
constant K [UJ is given by K = (7^(67/18 — C2) — 5n//9. Note that the color structure of T s 
involves only the factors CfCa and Cpnf^ all Cp terms cancel out. 




The two-loop soft anomalous dimension for e + e — > tt is plotted as a function of (3 with 
rif = 5 in Fig. 4. Note that Tg vanishes at /? = 0, the threshold limit, and diverges at (3 = 1, 
the massless limit. The inset in Fig. 4 shows more clearly the sharp increase of T s as j3 

(2) 

approaches 1. Figure 5 also displays T s but in a logarithmic plot which makes the small f3 
region more clear. 

It is instructive to study the small and large (3 behavior of T s in more detail. We first 
expand around (3 = and find 

rfL P = -^(3 2 {C F C A (18C 2 -A7) + 5C F n f ] + 0((3 i ). (7) 

(2) 

Note that T s is an even function of (3 and hence only even powers of (3 appear in the expansion. 
In Fig. 6 we see that just keeping the first term in the expansion as shown in Eq. ([7]) provides 
a good approximation to the complete result for small (3. The ratio of the expansion to the 

(2) 

complete T s remains within a few percent of unity for values of p up to around 0.3. Keeping 
more powers of (3 in the expansion results of course in a better approximation. Keeping powers 
through (3 12 results in an excellent approximation to over a wide range of f3, starting to 
differ by a few percent above (3 = 0.8 as seen from the corresponding curve in Fig. 6. 

As observed in Ref. [12], the two- loop soft anomalous dimension for processes with mass- 
less quarks is proportional to the one-loop result of Ref. [13] by the simple factor K/2, i.e. 
Fsmassiess = Tgrnassiess- This property was further explored and generalized in Refs. 

[T4"l [T5| [16] . We see this Kj 2 factor for heavy quarks as the first term in Eq. (J6j) . It is clear 
that for the massive case the situation is rather more complicated than for massless processes: 
although all the Cpiij terms of are included in (K/2) there are many additional CpC A 
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Figure 5: Logarithmic plot of T s for top quark production 
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terms. It is interesting to check however that in the massless limit of our T s this relation is 
satisfied. In Fig. 6 we also plot the ratio of (If/ 2) to r^. We see that indeed the ratio is 1 
at the limit (3 = 1, i.e. the massless limit, so the massless proportionality relation is recovered. 
However, for any other value of (3 the ratio is not one; in fact it is 1.144 near (3 = 0, the 
threshold limit, gradually decreases with increasing /3, and falls steeply to 1 in the region very 
close to (3 = 1. 

Given that the expansion around /?, Eq. (jTJ), gives very good approximations to 1^ at 
smaller (3 while the expression (K/2) Tg is a better approximation at large (3, we can derive 

(2) (2} 

an approximation to T s for all p values by starting with the p expansion of T s , and then 
adding {K/2) Fg and subtracting from it its (3 expansion: 



As seen from Fig. 6, adding just the (3 2 terms to (K/2) T s as in Eq. (jBJ) provides an excellent 

(2) 

approximation to the exact result for T s for all (3, never more than half a percent away from 
the exact value. It is quite remarkable that such a simple expression provides this good an 

(2) 

approximation to a rather long and complicated expression for T s . If we keep additional 
terms through order (3 12 then the approximation is truly impressive, not differing by more than 
one per mille anywhere in the (3 range, and in fact not differing by more than a few parts per 
million for most (3. At a linear collider with yfs = 500 GeV, (3 = 0.726 for a top mass of 172 
GeV. For this value of (3 the expansion through order (3 12 of Eq. (j7j) as well as the approximate 





1.05 



1.00 



0.95 



— 1 — 1 — 1 — 1 — 


— 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1- 


-1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


- 

- 


-- (K/2)r s (I) /r s (2) 
r <2) (0(p 2 ))/rf 

S approx x v " S 

— C pi Jo ( p 12 ))/rf 
C P (0(P 12 ))/r s (2) 

r s ( l P (°^ /r s (2) 


N 

N 

s 

\ 

\ 

\ - 
\ . 
t . 

1 






1 


1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 Ll—I 1 l_ 


\ 

\ 

\ 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



P 



Figure 6: Expansions and approximations to T s for e + e _ — > tt. 



expression at order (3 2 or higher of Eq. (jSJ) give excellent approximations. However, (K/2)T S ^ 
by itself is 11% larger than 

The result in Eq. (jf)]) can be written in terms of the cusp angle 7 = ln[(l +/?)/(! — f3)] as 
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where T^p = Cp{l coth 7 — 1), and is consistent with the cusp anomalous dimension of Ref. [5]. 

Finally we consider the two-loop soft anomalous dimensions for single top and for top pair 
(or bottom or charm pair) production at hadron colliders. Since the two-loop soft anomalous 
dimensions for massless quarks obey the simple proportionality relation to one loop, and the 
diagrams for single top and top pair processes involve eikonal graphs with both massless and 
massive eikonal lines, it is clear that for (3 near 1 the relation Tg = (K/2)Tg will hold for 
both single top and tt production. For smaller (3, however, this relation is no longer valid. 
For reference, at the Tevatron and the LHC the ti cross section receives most contributions in 
the region around 0.3 < (3 < 0.8 which peak roughly around (3 ~ 0.6. However, since even 
a 15% discrepancy between the exact Vg and (K/2) T s ^ may result in less than 1% effect 
in the total cross section it may turn out that using the massless proportionality relation is 
adequate for many purposes (see also the related discussion in [TT]). Nevertheless, a rigorous 
treatment requires that we include the extra terms for the massive case. More details will be 
given elsewhere. 
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